Force-Extension Relation and Plateau Modulus for Wormlike Chains 
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We derive the linear force-extension relation for a wormlike chain of arbitrary stiffness including 
entropy elasticity, bending and thermodynamic buckling. From this we infer the plateau modulus 
G° of an isotropic entangled solution of wormlike chains. The entanglement length L e is expressed 
in terms of the characteristic network parameters for three different scaling regimes in the entangled 
phase. The entanglement transition and the concentration dependence of G° are analyzed. Finally 
we compare our findings with experimental data. 
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Recently there has been increasing interest in biological 
material research Q. The physical properties of colloids, 
liquid crystals and macromolecular networks are of prime 
importance for the structure and function of biological 
entities such as cells and muscles. On the other hand 
biology provides physicists with some of the most perti- 
nent model systems to test their theories of soft matter 
M . Among these systems we will concentrate on macro- 
molecular networks here. These networks may be assem- 
blies of relatively flexible (DNA), semiflexible (actin) or 
rigid (microtubuli) molecules or complex compound sys- 
tems as in the case of the cytoskeleton of erythrocyte 
plasma membranes or in the glass body of the eye 
Especially the static and dynamic rheological properties 
of actin networks are crucial for an understanding of the 
mechanical stabilization and the motility of cells. 

Solutions of flexible polymers at high concentration are 
known to exhibit spectacular features in both their me- 
chanical response and molecular relaxation ffla], which 
are commonly attributed to the topological constraints 
due to the uncrossability of the polymers (entanglement). 
For time scales shorter than a characteristic time r the 
response of a solution of high molecular weight polymers 
to a periodic perturbation is elastic over an extended fre- 
quency interval ("rubber plateau") and resembles that 
of a permanently crosslinked network or gel. It is com- 
monly agreed that polymeric liquids form a temporary 
network, where entanglements play a similar role as per- 
manent crosslinks in gels. There have been several at- 
tempts to derive this successful phenomenological con- 
cept from other characteristic parameters of a polymer 
solution and thus gain some understanding of the un- 
derlying microscopic mechanisms HQ . Over the last few 
years experiments with actin [^|j9) revealed that the above 
qualitative picture also holds for semiflexible polymers. 
However, the theoretical understanding of entangled so- 
lutions of semiflexible polymers is much less developed 
than in the flexible case, and experimental data are often 
interpreted within the theoretical framework established 
for flexible coils or rigid rods respectively |9|,|iC 



In this letter we develop some basic concepts for a gen- 
eral theory of rheology of isotropic entangled solutions of 
wormlike chains. The wormlike chain |fl[ is the mini- 
mal model of an ideal (i.e. non self-avoiding) polymer 
of arbitrary stiffness. In contrast to the fractal Gaus- 
sian chain model 0] , which can serve only as an effective 
large scale model for rather flexible polymers, the worm- 
like chain model also faithfully reproduces the intrinsic 
stiffness of real polymers. The polymer is represented as 
a diffcrentiable space curve R s with its statistical prop- 
erties determined by the effective free energy 
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A central feature of this model is the inextensibility of 
the chain leading to the rigid constraint \dR s /ds\ = 1, 
which additionally has to be imposed on the contour. 
Due to the mathematical complications resulting from 
this constraint only few of the statistical properties of 
the wormlike chain can be extracted analytically, the 
most prominent being the mean square end-to-end dis- 



tance (R L ) 



— L 2 fj}(L/L p ), with the Debye function 



Jd{x) := 2{x — I + e~ x )/x 2 . (The persistence length L p 
is related to the bending modulus k by k = L p kgT.) For 
large L/L p this reduces to the power law (R L ) = 2L p L 
characteristic of a random walk of step length 2L p . Note 
that the model does not reproduce the swelling of real 
flexible coils, since the self- avoidance of the chain is ne- 
glected in the above effective free energy, Eq. ([!]). There- 
fore, the model is restricted to solutions with mesh size 
£ m not much larger than the persistence length L p of the 
polymers, in which case self-avoidance effects can safely 
be neglected. 

As a first step towards an understanding of the macro- 
scopic viscoelastic properties of an entangled network of 
wormlike chains one has to understand the elastic proper- 
ties of a single wormlike chain. The linear force-extension 
relation of a wormlike chain is obtained by the follow- 
ing argument. Consider a wormlike chain with one end 
clamped at fixed orientation at the origin. Apply a weak 
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force fn (directed along the unit vector n) to the other 
end [T^ ]. The configurational distribution function is 
then modified by a Boltzmann factor exp(/ni?i / k B T) . 
The extension SRl := ti((Rl) / — (Rl)) in the direction 
of the applied force to first order / is given by the linear 
extension coefficient fg '= d5R.L/df\f = o, 



fg^ = Jd s J ds' (cos 6 S cos 6 s i) 



ds (cos # s 



(2) 



By 6 S we denote the tilt angles of the tangents of the 
polymer contour with respect to n. The thermal average 
is to be taken under the constraint that the angle 9q at 
the clamped end is kept fixed. Standard methods |ll]] 
yield for f g "^ the dashed curves in Fig. |l|. In general, 
Eq. (0) predicts a polymer of contour length L to ap- 
pear more floppy if L p ~ L than in the high temperature 
limit (L p — > 0) and the low temperature limit (L p — > oo), 
when it contracts to a little ball or becomes a rigid rod 
respectively. In the flexible limit, where the chain be- 
comes an isotropic random coil, all curves fall together 
and reproduce entropy elasticity. But for stiff chains, as a 
consequence of the chain anisotropy, the force-extension 
relation depends strongly on the value of 6q. Obviously, 
9q = is an exceptional case. Whereas for all other an- 
gles 9q the ultimate asymptotic form of fg in the stiff 
limit is k/L 3 , at do = the force coefficient becomes 
/o — n 2 /k B TL 4 ; i.e. it is second order in the bending 
modulus and diverges at low temperatures T. The latter 
result was previously obtained in Ref. [^3) . Note that 9q 
is the angle between the applied force and the average 
orientation of Rl, i.e. for 6q = the force is parallel to 
Rl on average. Especially for T = the force is pulling 
or pressing on a rigid rod along its axis. In this limit the 
above expansion of the Boltzmann factor breaks down 
and we encounter the so called Euler buckling instability, 
i.e., the force-extension relation becomes highly nonlinear 
and the force coefficient in linear response does not exist. 
This situation is well known for foams and other cellular 
materials [[l4|. If we require JRl <C k B T, the buck- 
ling instability is evaded by thermal undulations, and 
we find a linear contribution to the force-extension re- 
lation ('thermodynamic buckling'). But with decreasing 
temperature the volume fraction k B TL/n (stored thermal 
energy over bending energy) occupied by the thermal un- 
dulations vanishes, and ultimately there remain no more 
undulations to be bent or pulled out, hence the diver- 
gence of /o with T^ 1 . We take as the force coefficient / 
of a 'general strand' of length L in a random network the 
average of fg 1 over all orientations 6q jig]: 



L 2 f ext (L/L p )/k B T, 



(3) 



with/ cxt (a;) := (2x — 3 + 4e _a: — e~ 2x )/3x 2 (the extension 
function). This result is also shown in Fig. [I]. In the stiff 
limit (L <C L p ) it reduces to / ~ k/L 3 . 
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. 1. The deformation SRl of a wormlike chain of given 
L to leading order in the applied force / as a function 
persistence length L p . See explanation in the text. 



Now, to relate the force-extension relation of the gen- 
eral strand to the observed elastic modulus of a polymer 
network in the rubber plateau regime we proceed in close 
analogy to solid state physics. In the harmonic approx- 
imation the elasticity tensor E of a monatomic Bravais 
lattice is written as E^m — —^2^ R -^RiDji(R)R k /2V, 
with V being the volume of the primitive cell, R the lat- 
tice vectors and D the matrix of second derivatives of 
the interaction potential with respect to lattice displace- 
ments. For an isotropic entangled polymer solution, we 
take the analogue of the primitive cell to be an entangle- 
ment volume V e and the analogue of the primitive vectors 
to be the average distance £ e between adjacent entangle- 
ments (in the embedding space). This scaling argument 
suggests that the storage modulus in the plateau regime 
should (up to a constant factor depending on the strain 
geometry) be given by G° ~ c e f e Q. Here f e is the force 
coefficient of polymer sections of length L e between ad- 
jacent entanglements and c e — V~ l their concentration. 
(As far as self-avoidance effects can be neglected, L e and 
£ e are related by the Debye function.) The situation may 
be visualized by external forces acting on contour ele- 
ments distributed with an average spacing L e along the 
polymer. Inserting / from Eq. (^ with L e substituted 
for L into the formula for G° we finally arrive at the fol- 
lowing explicit expression for the plateau modulus of an 
entangled solution of wormlike chains, 



G° ~ c e k B T 



.fp(L e /Lp) 

fext(L e /L p ) y C, 



c e k B T (L e > Lp) 
[L e < Lp) . 



(4) 



The entanglement length L e is obviously the crucial 
quantity in Eq. (||). In the literature several scaling ideas 
H have been reported on how L e and £ e may be derived 
from the known static properties of a flexible polymer 
network. Note however that Eq. (^) holds independently 
of such considerations. For a homogeneously crosslinkcd 
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gel of semiflexible or rod-like polymers £ e can essentially 
be identified with the mesh size £ m of the network. In 
this case Eq. (@) predicts G° oc kc 2 in the stiff limit. 
We conjecture that for a solution of wormlike chains of 
arbitrary stiffness one has to distinguish three different 
regimes. We will treat the limiting cases of scale invariant 
chain structure - i.e. a virtually Gaussian or straight con- 
formation respectively - in a very similar manner. The 
breaking of scale invariance due to Eq. (0) gives rise to 
an intermediate regime for chains with L ss L p , which 
will be discussed subsequently. 

For a weakly bending contour we have from Eq. (|l|) 
the scaling relation R^ 2 = 2L 3 /3L p for the transverse 
amplitudes of the largest bending undulations. If 
these amplitudes are smaller than the mesh size £ m of 
the surrounding polymer network, i.e. £^ > 2L 3 /3L p , 
then the bending undulations are not substantially per- 
turbed and are supposed to be rather irrelevant to the 
question of entanglement. In this case we should thus be 
allowed to represent the polymers as straight (but not 
rigid) "rods" in our derivation of L e . In the opposite 
extreme [L,^ m 3> L p ) of a strongly coiled polymer con- 
formation the polymer may be represented by a fractal 
curve (or a freely jointed chain of "blobs" , if screened self- 
avoidance is to be included) . We feel that the flexible case 
has been described successfully earlier and will adapt 
this approach to straight rods now. It is based on the cru- 
cial observation that polymer ends are not contributing 
efficiently to long-lived entanglements. Namely, if entan- 
glements would depend on dangling ends, they could not 
be long-lived as compared to unperturbed, free fluctua- 
tions of the polymer and hence could not give rise to a 
rubber plateau. To be specific, we assume that an entan- 
glement requires a sufficient number of non-end neighbor- 
ing polymer segments that on average restrict the lateral 
degrees of freedom of a test chain. Consider a sphere of 
radius £ e around such a mean entanglement point. Then 
for a given monomer concentration c and volume fraction 
cv the excluded volume in this "primitive cell of entan- 
glement" of volume V e = 47r(£ e /2) 3 /3 is given by cvV e . 
In order to achieve an entanglement one requires that a 
certain amount of polymer material, Cir{a/2) 2 L e , is con- 
tained in the test volume. Here a denotes the lateral 
diameter of the polymer. The quantity C is a geometry 
factor which measures the amount of polymer material in 
the test volume, if K polymers cross the sphere around 
the test chain [|l8| . We determine L e by equating the ex- 
cluded volume (reduced by the contribution coming from 
free ends) with the volume of polymer material needed 
for an entanglement 
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The implicit equation Eq. (^) can be solved analytically 
in the random coil limit |H| and for a straight conforma- 
tion. For the latter we find 



where := a^/3C/2cv ~ VC£ m is the entanglement 
length in the limit of infinitely long molecules. Eq. 
(g) describes an entanglement transition of the poly- 
mer solution characterized by a cusp singularity of the 
entanglement length L e as a function of the polymer 
length L or the volume fraction cv. The phase bound- 
ary between the entangled and the disentangled regime is 
given by either of the two equations c*(L) = 27Cc(L)/4, 
L*(c) = 3v3L£°(c)/2, where c denotes the geometrical 
overlap concentration cv = 3a 2 /2L 2 . The value of the 
entanglement length and the contour length at the cusp 
singularity are related by L* = 2L*/3. The above re- 
sults have some important consequences on the rheolog- 
ical properties of semiflexible polymer solutions. Upon 
taking the experimental value for the geometric constant 
C = 9.1 Q of flexible polymers to be a universal quantity 
(also valid in the rod-limit) one estimates that the critical 
polymer length for the solution to show entanglement has 
to be about eight times the mesh size. Accordingly, the 
critical concentration c* is predicted to be almost two 
orders of magnitude larger than the overlap concentra- 
tion c, c*/c = 27C/4. In the intermediate concentration 
regime, c* > c > c, there is already a significant overlap 
of the semiflexible polymers but no long-lived entangle- 
ments leading to a rubber plateau regime. In this disen- 
tangled phase the magnitude of the storage modulus is 
supposed to show a linear concentration dependence. 

Finally we comment shortly on the intermediate case of 
a network of wormlike chains with a mesh size £ m smaller 
than the persistence length L p and the amplitudes of the 
largest bending undulations i?jh To distinguish it from 
the case discussed above, which could be called the "rod- 
like" regime, we will address it as the "snake-like" regime. 
It is characterized by the property that all bending un- 
dulations with wavelength longer than a critical wave- 
length, the "deflection length" jl7) A ~ {ZL^/2) 1 / 5 , 
are perturbed by the network. For the snake-like regime 
we thus identify the entanglement length for an in- 
finitely long polymer with A. We expect the qualitative 
features of the entanglement transition derived above for 
the rod-like regime to hold also in the snake-like regime. 
The isotropic entanglement volume V e ~ £g in Eq. (||) 
has now to be replaced by V e ~ £ e i?£ 2 . The implicit 
dependence of V e on L p again reflects the broken scale 
invariance in the snake-like regime. We leave the prob- 
lem of the crossover between the snake-like and the scale 
invariant cases for further investigation. 

Now we turn to the comparison of our results with 
available experimental data. We suppose that the ex- 
istence of a disentangled phase above the overlap con- 
centration c, as predicted by Eq. (|B|), is likely to explain 
some discrepancies of the Doi-Edwards theory |||l9| for 
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the rotational diffusion of rigid rods in a semidilute so- 
lution with experimental data (20). The Doi-Edwards 
theory is based on the assumption that for concentra- 
tions larger than the overlap concentration c there is a 
separation of time scales. One anticipates that each step 
of the rotational diffusion process is determined by the 
constraint that each rod is confined to remain within an 
angular range /L during the time it takes a rod to dif- 
fuse a distance equal to its length. The tube of radius 
£t oc 1/cL is assumed to impose a long-lived topologi- 
cal restriction on the motion of the rods. However, as 
we have argued above, long-lived entanglements emerge 
only at a much higher concentration c* = 27Cc/4. This 
would explain why the onset of entanglement as defined 
by a marked decrease in the rotational diffusion coeffi- 
cient occurs at a concentration c CX p fully two orders of 
magnitude above the overlap concentration |2(i(| . 

As an important practical application of the above 
ideas we already mentioned actin, a semiflexible macro- 
molecule, which is of major biological interest but is also 
an almost ideal model system for physicists 21 1 . It is well 
suited to test our ideas, because it is characterized by a 
large ratio L p /a (~ 10 3 ) and thus by a broad semidilute 
regime, so that the wormlike chain model applies with- 
out modification over several orders of magnitude in con- 
centration. The average length of the molecules can be 
adjusted by adding so called actin binding proteins such 
as gelsolin or severin. Existing data on actin |p|,p|,p2[ give 
only an incomplete picture of the rather complex situ- 
ation sketched above but seem to support our results. 
For short (rod-like) filaments a length dependence of the 
plateau modulus near the entanglement transition has 
been observed [|j , which is qualitatively well described by 
Eq. (Q) and (^), but is somewhat smeared out (probably 
as an effect of sample polydispersity). In the entangled 
phase Eq. (fij) together with Eq. (||) predicts a plateau 



modulus G° 



oc KC 



for the rod-like case far from the en- 



tanglement transition. Near the transition the concentra- 
tion dependence is enhanced according to Eq. (^) . In the 
snake- like case as defined above we have L e ~ (£f)£m) 
with £ TO cx c" 1 / 2 and hence G° cx c 5 / 3 ■K 1 / 3 (k B T) 2 / 3 . 
In experiments with actin the exponent of the observed 
power law for G°(c) ranges from 1.7 to 2.3 in the entan- 
gled regime |§|, |l3| , p2f . As a critical test of our ideas we 
suggest a comparison of the plateau modulus for actin 
networks with and without tropomyosin, which is known 
to cause a considerable stiffening of actin filaments. The 
rod-like regime and the snake-like regime should be read- 
ily discernible due to their markedly different dependence 
of G° on k. 

In summary, we have derived the force-extension rela- 
tion for a wormlike chain and discussed some of its con- 
sequences for the viscoelastic properties of entangled so- 
lutions and gels of semiflexible polymers. Especially, we 
analyzed the entanglement transition and predicted vari- 
ous exponents for the dependence of the plateau modulus 



G° on concentration and bending rigidity. 
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